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DESCRIPTION 



X-RAY TOMOG RAPH 



TECHNICAL FIELD 



The present invention relates to a tomograph 



which generates a tomographic image of an examinee 
using projection data obtained from a radiation source 
5 moving in the body axis direction relative to the 
examinee through a radiation detector. 

BACKGROUND ART 



projection method will be explained. A Feldkamp 
10 method, Wang method, IHCB method and Pi-method proposed 
as three-dimensional back projection methods are three- 
dimensional back projection methods which capture a 
cone beam spreading (having an angle of inclination) in 
both the slice (body axis) direction and channel 
15 (rotation)' direction irradiated to multi-row radiation 
detectors as a collection of a plurality of rows of fan 
beams spreading only in the channel direction, carry 
out filter correction processing similar to a two- 
dimensional back projection method on the fan beam 
20 projection data obtained from each detector row or 

parallel beam projection data obtained by replacing the 
fan beam by parallel beam through rearrangement 
processing and carry out back projection processing 



A conventional three-dimensional back 



along the trace of the beam to thereby obtain a 
reconfigured image . 

Fig. 7 shows a reconf igurable condition of a 
Wang method and Fig. 8 shows a reconf igurable condition 
of a Pi-method. Here, reference character FOV denotes 
an effective field of view region, SOD denotes a 
distance between an X-ray tube and the go-around axis 
of a CT device and SID denotes a distance between the 
X-ray tube and CT device. The Wang method is a method 
corresponding to a Feldkamp method adapted to image 
taking of a spiral orbit and has a back projection 
phase width of n to 2n. 

An example of a Pi-method is disclosed in JP- 
A-ll-253434. This is also a back projection method 
applicable to image taking of a spiral orbit and is a 
reconfiguration method for back projecting a n range in 
which the phase varies from one voxel to another to 
improve a bed moving speed using the Wang method. The 
Pi-method can set the back projection phase range for 
each voxel to n by limiting the vertical direction of 
an X-ray beam to be back projected using a spiral 
located opposite to the X-ray focal position. 

An example of the IHCB method is disclosed in 
JP-A-11-4823 . This method consists of an algorithm for 
back projecting a back projection phase range which 
varies from one voxel to another and the back 
projection phase width is either n or an entire 
possible data range which varies from one voxel to 



another . 

Next, problems of these conventional 
technologies will be explained. 

The Feldkamp method is an image 
reconfiguration method for image taking of a circular 
orbit and is not applicable to image taking of a spiral 
orbit. The Wang method is an image reconfiguration 
method for image taking of a spiral orbit and can 
correct influences of movement of an examinee, which is 
practiced by the conventional two-dimensional back 
projection method by extending the back projection 
phase width beyond n (increasing data redundancy) , but 
results in a poor data utilization rate and the pitch 
(hereinafter referred to as "measuring throughput" ) of 
the spiral during image taking needs to be very small. 
By improving the Pi-method and IHCB method so that the 
back projection phase range according to the Wang 
method is widened, their respective measuring 
throughputs can be drastically improved compared to the 
Wang method, but they are the back projection methods 
within the n range with data redundancy completely 
eliminated, and therefore data may be discontinuous at 
the start phase and end phase of the back projection 
phase range due to influences of movement of the 
examinee, which is likely to become a strong artifact 
and appear on the image. 

Here, data redundancy will be explained. The 
data redundancy refers to a breadth of a phase range 



within which not only phase data but also opposed pha 
data is acquired. According to a three-dimensional 
back projection method, data redundancy changes from 
one voxel to another. For example, as shown in Fig. 
22, when back projection is performed from data 
obtained by rotating the phase of a radiation source 
180 degrees, the contributing data phase range change 
from one reconfiguration pixel to another and a pixel 
has data having a phase range of 180 degrees or more, 
while a pixel b can only acquire data of 180 degrees < 
less. Furthermore, it is also necessary to consider 
the beam width in the body axis direction and in this 
way data redundancy changes from one pixel to another 
in a complicated manner. For this reason, a 
complicated redundancy correction is required. 

One of problems of these conventional three- 
dimensional reconfigurations is an increase in a 
calculation time. 

Therefore, when an increase in the amount of 
calculation from a parallel beam two-dimensional back 
projection method to a parallel beam three-dimensional 
back projection method is analyzed, the increased 
calculation causes (1) an increase in the number of 
times one-dimensional rearrangement processing is 
performed, (2) an increase in the number of times 
reconfiguration filter processing is performed and (3) 
an addition of calculation of detector row addresses 
during back projection processing. Here, the main 



processing that occupies the calculation time in the 
two-dimensional back projection method and three- 
dimensional back projection method is back projection 
processing. 

The loads of calculation of the distance 
between the focus and reconfiguration point during the 
calculation of detector row addresses and arcsin 
calculation (calculation of the z position of the focus 
of the parallel beam of the following Expression 1) are 
particularly large and occupy the major portion of 
causes of increases in the calculation time. 
[Expression 1] 

z S = (J * (0 + arcs in (t^SOD) ) / 2 n) + z so 
See Fig. 29. 

Suppose SOD is a distance between a radiation 
source and a go-around axis, <|> is a phase angle of the 
parallel beam, J is a relative movement distance from a 
radiation source to an examinee per rotation of a 
scanner on a radiation detector 13, ti is the position 
in the channel direction, z s is the position of the 
radiation source 11 in the z direction and z S o is z s 
when the go-around phase of the radiation source is 0. 
Therefore, if these calculations can be simplified, it 
is possible to significantly increase the processing 
speed of the tomograph. 

It is an object of the present invention to 



provide a tomograph capable of suppressing generation 
of the distortion attributed to data discontinuity and 
obtaining a tomographic image of high image quality not 
eliminating data redundancy but rather using it in 
three-dimensional back projection calculations. 

It is another object of the present invention 
to provide a tomograph capable of simplifying arcsin 
calculation on a fan-parallel beam conversion and back 
projection processing according to a set FOV range in 
three-dimensional back projection calculations and 
significantly increasing the processing speed of the 
tomograph without degrading image quality. 

DISCLOSURE OF THE INVENTION 

1. In order to attain the above described 
objects, the present invention is a tomograph 
comprising a radiation source and a radiation detector 
arranged opposite to each other, between which a bed 
with an examinee placed thereon is provided, the 
radiation source and radiation detector turning around 
the bed which can be moved with respect to this go- 
around axis, radiation irradiated from the radiation 
source and passing through the examinee being detected 
using the radiation detector, and reconfiguration means 
for creating a three-dimensional tomographic image in a 
region in concern of an object from the detected 
projection data, wherein the reconfiguration means 
determines for each voxel a projection data range 



capable of back projection having an operating 
projection data phase width of 180 degrees or more, 
superimposes a reconfiguration filter, assigns weights 
to data of the same phase or opposite phase for each 
phase for this projection data range and performs 
three-dimensional back projection on this filter- 
processed projection data over the determined data 
range capable of back projection along the irradiation 
trace of the radiation beam. 

Since the tomograph of the present invention 
determines the projection data phase range used for 
each voxel, it is possible to determine the projection 
data phase range for each voxel so that absolute values 
of the angles of inclination of radiation beams become 
the same at both ends of the projection data, thereby 
use projection data with a small cone angle, provide 
redundancy using weighting means and correct the data 
for each voxel using a weighting function, thereby 
suppressing generation of the distortion attributed to 
discontinuity in the data phase direction and obtain 
images of high quality. The tomograph of the present 
invention requires no redundancy processing which would 
require complicated calculations, thus making it 
possible to create images at high speed. 

2. The present invention described in the 
item 1 is characterized in that when determining the 
above described data range, a projection data range is 
determined so that the difference in the absolute 
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values of cone angles at both ends of the projection 
data range used is reduced. 

3. The present invention described in the 
item 2 is characterized in that the projection data 
phase width used is determined so as to be the same 
phase width for each voxel. 

The tomograph according to the invention 
described in the items 2 and 3 is characterized in that 
determining means for determining the projection data 
phase range used for each voxel determines the 
projection data range so that the difference in the 
absolute value of cone angles at both ends of the 
actually used projection data range becomes small or 
determines the projection data range so that the 
projection data phase width used has the same phase 
width for each voxel, which allows projection data with 
a small cone angle to be used. Furthermore, by 
equalizing the absolute values of angles of inclination 
of radiation beams at both ends of the projection data 
exactly, it is possible to calculate the position of 
the detector row direction from the data start 
direction or end direction simultaneously and further 
calculate the same phase range at the time of back 
projection of each reconfigured voxel and thereby 
determine a weighting function for redundancy 
corrections using a single expression and perform 
calculations at high speed. 

4. The present invention described in the 



item 1 is characterized in that the projection data 
range capable of back projection is either 270 degrees 
or 360 degrees . 

The tomograph of the present invention 
described in the item 3 uses either 270 degrees or 360 
degrees as the projection data range capable of back 
projection, and assigns weights to data using 270 
degrees in the phase direction, and can thereby reduce 
discontinuity at the data end to a minimum. This 270- 
degree data corrects a discontinuity at the 180-degree 
data end using a data phase with smallest discontinuity 
having a 90-degree phase difference and can reduce data 
discontinuity to a minimum, and thereby realize 
reconfiguration of high quality. 

5. The invention described in any one of the 
items 1 to 4 is characterized in that projection data 
whose number of images taken per rotation is a multiple 
of the number of sides C of a polygonal display pixel 
is acquired and the reconfiguration means comprises 
back projection means for superimposing the 
reconfiguration filter on this projection data, 
grouping data at the same channel position and having 
projection phases in the go-around direction shifting 

by Ntt/2 (N=l, 2, 3, ) [rad] at a time and 

performing back projection to a square image array 
group by group. 

6. The invention described in any one of the 
items 1 to 4 is characterized in that the 



reconfiguration means converts the projection data 
obtained to data including fan beam data and parallel 
beam data whose number of images taken per rotation is 
a multiple of the number of sides C of a polygonal 
display pixel, superimposes the filter on this 
projection data, groups data at the same channel 
position and having projection phases in the go-around 

direction shifting by Ntc/2 (N=l, 2, 3, ) [rad] 

at a time and performs back projection to a square 
image array group by group. 

The tomograph of the invention described in 
the items 5 and 6 is a method for enhancing the speed 
of back projection requiring the maximum calculation 
time in creating an image. In order to enhance the 
speed of back projection, the present invention takes 
advantage that the shape of the reconfigured image 
array is polygonal and that image taking is performed 
while circling around the reconfigured image, the 
invention described in the item 5 takes images with a 
view which is a multiple of the number of sides of a 
display pixel, performs fan beam reconfiguration and 
the invention described in the item 6 converts data to 
data whose number of views is a multiple of the number 
of sides of a display pixel through rearrangement 
processing and performs parallel beam reconfiguration. 
In all cases, the invention groups projection data 
whose phase in the go-around direction shifts by 
Ntc/2 (N=l, 2, 3, ) [rad] at a time, back projects 



the square image group by group, and can thereby reduce 
the number of times the channel direction position in a 
full reconfiguration and interpolation coefficient are 
calculated. This is because when the reconfigured 
image is square, the data of a phase differing exactly 

by Ntc/2 (N=l, 2, 3, ) [rad] and the square 

reconfigured image have the same positional 
relationship. Furthermore, the number of views is set 
to a multiple of 4 to calculate data of a phase 

differing by Ntc/2 (N=l, 2, 3, ) [rad] exactly and 

it is possible to create images by calculating channel 
positions within a range of 1/4 of a full revolution 
(tc/2 [rad] ) in the cases of both a full reconfiguration 
and a half reconfiguration. In this way, in the case 
of a full reconfiguration, the amount of calculation 
becomes 1/4 and though calculations are carried out 
using one calculator, a calculation result close to a 
result of a parallel calculation using four calculators 
can be obtained and it is possible to realize high 
performance at low cost. 

7. The invention described in any one of the 
items 1 to 6 is characterized in that associating means 
is provided for associating pixel intervals in the body 
axis direction of the image using polygonal display 
pixels with the relative moving speed between the 
object and the radiation source in the go-around axis 
direction . 

8. Furthermore, the invention described in 



the item 7 is characterized in that the associating 
means is constructed so that the relationship between 
pixel interval rpitch in the body axis direction of the 
square image and the relative moving speed J in the go- 
around axis direction of the object and the radiation 

source is expressed by J=2-N-rpitch (N=l, 2, 3 ). 

9. The tomograph according to the invention 
described in the items 7 and 8 is characterized in that 

at the phase of Ntt (N=l, 2, 3, ) [rad] of the 

radiation source, the position on the radiation 
detector at which the beam passing through a voxel I 
(x, y, z) whose body axis direction position is Z [mm] 
and a voxel I (-x, -y, NJ/2+Z) whose body axis 
direction position is- N • J/2 + Z [mm] intersects remains 
the same, and therefore when a beam passing through a 
voxel is calculated at a certain view at the time of 
back projecting, this is equivalent to simultaneous 
calculations of the row positions of phases differing 

by Ntc (N=l, 2, 3, ) [rad] from each other and 

when an image is generated from data with a plurality 
of revolutions obtained by taking images through spiral 
scanning, it is possible to enhance the speed of back 
projection which requires a maximum time for image 
generation . 

10. Furthermore, the present invention is a 
tomograph comprising a radiation source and a radiation 
detector made up of two-dimensionally arranged 
detection elements, arranged opposite to each other, 



between which a bed with an examinee placed thereon is 
provided, the radiation source and radiation detector 
turning around the bed which can be moved with respect 
to this go-around axis, radiation irradiated from the 
radiation source and passing through the examinee being 
detected using the radiation detector, and 
reconfiguration means for creating a three-dimensional 
tomographic image in a region in concern of the 
examinee from the detected projection data, wherein the 
reconfiguration means determines a projection data 
phase range capable of back projection for each 
reconfigured voxel, calculates an approximate straight 
line for a curve indicating the radiation source 
position with respect to the channel direction position 
of parallel beam projection data corresponding to the 
region in concern obtained by a parallel beam of a 
parallel shape viewed from the go-around axis direction 
generated from the radiation source, corrects each row 
of the projection data by multiplying a coefficient 
which is dependent on the angle of inclination of 
radiation from the radiation source, carries out one- 
dimensional rearrangement processing for obtaining 
parallel beam projection data from the fan beam 
projection data obtained from a fan-shaped fan beam 
viewed from the go-around axis direction generated from 
the radiation source, and superimposes the 
reconfiguration filter on the parallel projection data 
to generate filter-processed parallel projection data 



and three-dimension back projects the parallel beam 
projection data subjected to the filter processing 
based on the determined projection data range capable 
of. back projection to the back projection region 
corresponding to the region in concern along the 
approximate irradiation trace using the approximate 
straight line . 

The tomograph of the present invention three- 
dimension back projects the filter-processed parallel 
beam projection data to the back projection region 
corresponding to the region in concern based on the 
projection data range capable of back projection 
determined by the operating data phase range 
calculating means along the approximate irradiation 
trace of the radiation beam calculated using an 
approximate straight line by the approximate straight 
line calculating means for calculating an approximate 
straight line for the curve indicating the radiation 
source position relative to the channel direction 
position corresponding to the region in concern of the 
parallel beam projection data obtained by the parallel 
beam of a parallel shape viewed from the go-around axis 
direction generated from the radiation source, and 
therefore as opposed to the conventional focus position 
calculation of a parallel beam which includes arcsin 
calculation and has an increased load, this arcsin 
calculation is replaced by an approximate straight 
line, simplifying the amount of calculation in the 



parallel beam three-dimensional back projection method 
and thereby significantly increasing the processing 
speed of the tomograph. 

11. The present invention described in the 
item 10 is characterized in that the reconfiguration 
means performs redundancy correction weighting for 
generating a weighting factor from a weighting function 
in the phase direction to correct data redundancy at 
each phase according to the phase width of this 
determined projection data, and the parallel beam 
three-dimensional back projection means assigns the 
weighting factor obtained by the redundancy correction 
weighting means to the projection data within the 
determined projection data phase range and performs 
three-dimensional back projection along the approximate 
trace to the back projection region. 

12. The tomograph of the present invention 
described in the item 11 is characterized in that in 
determining the projection data phase range, it is 
possible to determine the phase range of fit [rad] in 
the view direction and perform a redundancy correction 
using the weighting function by the redundancy 
correction weighting means. Thus, it is possible to 
provide data with redundancy (extending the back 
projection phase width beyond 180 degrees), assign 
weights using the weighting function, reduce 
discontinuity at the data ends (at the start/end of 
image taking) and obtain an image with the influence of 
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movement of the examinee reduced to a minimum. 

13. Furthermore, the invention described in 
the item 10 is characterized in that the operating data 
phase range calculation means determines the projection 
data range capable of back projection for each 
reconfigured voxel so that the maximum cone angle of 
the beam back projected for each voxel becomes 
narrowest . 

The tomograph of the invention described in 
the item 13 determines the back projection phase range 
for each voxel by the operating data phase range 
calculation means so that the maximum cone angle 
becomes a minimum, and can thereby reduce influences of 
deterioration of image quality by the cone angle to 
obtain better image quality and improve the relative 
moving speed (so-called measuring throughput) between 
the examinee and focus in the Z direction. 

14. Furthermore, the invention described in 
the item 10 is characterized in that in calculating the 
operating data phase range calculation, the projection 
data range capable of back projection for each 
reconfigured voxel is determined so that the phase 
direction range of the beam back projected for each 
voxel is set to the narrowest possible range. 

The tomograph of the invention described in 
the item 14 determines the back projection phase range 
for each voxel by the operating data phase range 
calculation means so that the number of views becomes 



small, and can thereby improve time resolution for each 
voxel. Furthermore, combining the invention with the 
redundancy correction weighting means described in the 
item 13 can obtain better image quality in a region 
where the examinee moves fast. Furthermore, by setting 
the back projection phase range for each voxel to the 
time range in which images are taken at the same time 
whenever possible so that the time position of the 
respective voxels in the displayed images come closer 
to one another, it is possible to shorten the time 
width contributing to the reconfigured image and 
improve time resolution. 

BRIEF DESCRIPTION OF DRAWINGS 

Fig. 1 is a schematic view of a tomograph of 

the present inventions- 
Fig. 2 is a block diagram of the tomograph 

shown in Fig. 1; 

Figs. 3A and 3B are conceptual diagrams 

showing a focus trace of a circular orbit scan and 

spiral orbit scan; 

Figs. 4A and 4B are side views of the waist 

part of a single radiation detector and multi-row 

radiation detector; 

Figs. 5A and 5B illustrate a collimation 

thickness of an X-ray beam per row of the single 

radiation detector and multi-row radiation detector; 

Fig. 6 is a flow chart showing a processing 
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operation by general reconfiguration means; 

Fig. 7 is a plan view showing a positional 
relationship between the radiation detector 13 and beam 
as a reconf igurable condition according to a Wang 
method; 

Fig. 8 is a plan view showing a positional 
relationship between the radiation detector 13 and beam 
as a reconf igurable condition according to a Pi-method; 

Fig. 9 is a flow chart showing a processing 
operation of an embodiment of the reconfiguration means 
of the present invention; 

Fig. 10 is a flow chart showing an operation 
of the operating data position range determining 
processing shown in Fig. 9; 

Figs. 11A and 11B are a perspective view and 
an exploded view showing spiral traces of the radiation 
source and radiation detector; 

Fig. 12A and 12B are a perspective view and 
an exploded view illustrating the operation of the 
rearrangement processing shown in Fig. 9; 

Figs. 13A and 13B are a perspective view and 
an exploded view illustrating other rearrangement 
processing; 

Fig. 14A is a spiral measuring diagram when 
180-degree data is used; 

Fig. 14B is a characteristic diagram showing 
a weighting function corresponding to spiral 
measurement when 180-degree data is used; 



Fig. 15A is a spiral measuring diagram when 
270-degree data is used; 

Fig. 15B is a characteristic diagram showing 
a weighting function corresponding to other spiral 
measurement when 270-degree data is used; 

Fig. 16 is a flow chart showing a processing 
operation of another embodiment of the reconfiguration 
means of the present invention; 

Fig. 17 is a perspective view showing an 
example of a back projection data phase range when 180- 
degree phase range (f=l) is used; 

Fig. 18 is a perspective view showing an 
example of a back projection data phase range when a 
phase range (Kf<2) from 180 degrees to 360 degrees is 
used; 

Fig. 19 is a perspective view showing an 
example of back projection data phase range when 360- 
degree phase range (f=2) is used; 

Figs. 20A and 20B are characteristic diagrams 
of a weighting function illustrating the redundancy 
weighting processing shown in Fig. 16; 

Figs. 21A to 21C are characteristic diagrams 
of a weighting function for each phase illustrating 
redundancy weighting processing shown in Fig. 16; 

Fig. 22 is a plan view showing a phase range 
capable of back projection; 

Fig. 23 is a flow chart showing a processing 
operation of another embodiment of the reconfiguration 



means of the present invention; 

Fig. 24 is a flow chart showing a processing 
operation of another embodiment of the reconfiguration 
means of the present invention; 

Fig. 25 illustrates back projection in group 
units in the grouping processing shown in Fig. 23; 

Fig. 26 illustrates another back projection 
in group units in the grouping processing shown in Fig. 
23; 

Fig. 27 is a flow chart showing a processing 
operation of a further embodiment of the 
reconfiguration means of the present invention; 

Fig. 28A, 28B and 28C show a relationship 
between the figure illustrating calculation processing 
of the approximate straight line shown in Fig. 27, cone 
angle, X-ray source and reconfiguration; 

Fig. 29 illustrates determining processing of 
the projection data phase range; 

Fig. 30 is a flow chart showing a processing 
operation of a still further embodiment of the 
reconfiguration means of the present invention; 

Fig. 31 is a flow chart showing an operation 
of phase range calculation processing on the data used 
in Fig. 30; and 

Fig. 32A to Fig. 32D show a relationship 
between a fan beam and parallel beam. 



BEST MODE FOR CARRYING OUT THE INVENTION 



Embodiment 1 

With reference now to the attached drawings, 
embodiments of the present invention will be explained 
in detail below. 

Fig. 1 is a schematic outside view of a 
tomograph according to an embodiment of the present 
invention. The tomograph comprises a scanner 1 used 
for image taking, a bed 2 to place and move an 
examinee, an input device 3 made up of a mouse and 
keyboard, etc., for inputting measuring reconfiguration 
parameters such as bed moving speed information and 
reconfiguration position, a calculation device 4 for 
processing data obtained from a multi-row detector and 
a display device 5 which displays a reconfigured image. 

Fig. 2 is a block diagram of the tomograph 
shown in Fig. 1. 

The scanner 1 consists of a bed 2, a high 
voltage switching unit 8, a high voltage generation 
device 9, a radiation source 11 such as a radiation 
generation device having a radiation control device 10, 
a radiation detector 13 placed opposite to the 
radiation source 11 with respect to an examinee 12, a 
go-around drive device 14 which drives this radiation 
detector 13 and radiation source 11 in the go-around 
direction and a collimator 15 which controls a 
radiation region to be irradiated from the radiation 
source 11, etc. The scanner 1 further consists of a 
collimator control device 16 which controls the 
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collimator 15, a scanner control device 17 which 
controls the go-around drive device 14, a bed movement 
measuring device 19 which measures an amount of 
relative movement between the bed control device 18 
which controls the bed 2 and bed 2 and a central 
control device 20 which controls these devices. 

Image taking conditions (bed moving speed, 
tube current, tube voltage, slice position, etc.) and 
reconfiguration parameters (region in concern, 
reconfigured image size, back projection phase width, 
reconfiguration filter function, etc.) are input from 
the input device 3, a control signal necessary for 
image taking is sent from the central control device 20 
to the radiation control device 10, bed control device 
18 and scanner control device 17 based on the 
instruction and upon reception of an image taking start 
signal, image taking is started. When image taking is 
started, the radiation control device 10 sends a 
control signal to the high voltage generation device 9, 
a high voltage is applied to the radiation source 11, 
and radiation is irradiated from this radiation source 
11 to the object 12. At the same time, a control 
signal is sent from the scanner control device 17 to 
the go-around drive device 14, and the radiation source 
11, radiation detector 13 and preamplifier 21 turn 
relative to the object 12. On the other hand, the bed 
2 carrying the examinee 12 is stopped by the bed 
control device 18 during a circular orbit scan or 



translated in parallel in the go-around axis direction 
of the radiation source 11, etc., during a spiral orbit 
scan. The go-around drive device 14, scanner control 
device 17 and bed control device 18, etc., constitute a 
drive device which turn the radiation source 11 and 
radiation detector 13 relative to the- examinee 12 and 
which is relatively movable in the axial direction of 
the examinee 12. 

With the irradiation region restricted by the 
collimator 15, the radiation irradiated from the 
radiation source 11 is absorbed and attenuated by each 
tissue inside the examinee 12, passed through the 
examinee 12 and detected by the radiation detector 13. 
The radiation detected by the radiation detector 13 is 
converted to a current, amplified by the preamplifier 
21 and input to the calculation device 4 as a 
projection data signal. The projection data signal 
input to the calculation device 4 is processed by the 
reconfiguration means 22 for reconfiguring the image 
inside the calculation device 4. 

The reconfigured image is saved in a storage 
device 23 in an I/O device 50 and displayed by an image 
processing device 26 on the display device 5 as a 
tomographic image . 

Figs. 3A and 3B are conceptual diagrams 
showing focus orbits of a circular orbit scan and 
spiral orbit scan. 

Fig. 3A shows a movement trace 24a of the 



radiation source (focus) during a circular orbit scan 
and Fig. 3B shows a movement trace 24b of the radiation 
source (focus) during a spiral orbit scan. If the 
detector is formed of a single row, when images are 
taken on a circular orbit as in the case of the 
movement trace 24b, it is possible to accurately 
reproduce the images at the positions of the radiation 
source by carrying out filter correction two- 
dimensional back projection. However, when images are 
taken on a spiral orbit as in the case of the movement 
trace 24b, carrying out filter correction two- 
dimensional back projection alone results in streak- 
shaped artifact at that position due to data 
discontinuity at the end position of image taking. 
Thus, by applying data interpolation to the data 
obtained on the spiral orbit as in the case of the 
movement trace 24b, the data is corrected to the 
circular orbit data like the movement trace 24a and 
then filter correction two-dimensional back projection 
is carried out. In this way, it is possible to obtain 
an image with reduced discontinuity. The degree of 
artifact in this case is determined by the degree of 
discontinuity in the X-ray source trace, that is, the 
degree of artifact changes depending on the moving 
speed of the examinee. For example, in a single row 
type spiral scanning X-ray tomograph (SDCT), the spiral 
pitch (ratio of the moving speed of the examinee to the 
thickness of the X-ray beam in the go-around axis 



direction) is generally used to an extent that 
substantially the entire image taking region can be 
covered with the data on the opposite side taken into 
consideration . 

Figs. 4A and 4B are schematic side views of 
the single row radiation detector 13a and multi-row 
radiation detector 13b. 

In Fig. 4B, a plurality of multi-row 
radiation detectors 13b whose width per row is narrower 
than that of the single row radiation detector 13a in 
Fig. 4A are arranged in a plurality of rows in the go- 
around axis direction, realizing a wider detector than 
the single row radiation detector 13a as a whole. 

Figs. 5A and 5B are schematic side views 
illustrating a thickness in the go-around axis 
direction (hereinafter, referred to as a "detector 
collimation thickness DCT") of the radiation beam per 
one row of detectors at the position of the collimator 
15 in the case of using the single row radiation 
detector 13a and the multi-row radiation detector 13b, 
respectively . 

In the case of the multi-row radiation 
detector 13b shown in Fig. 5B, the detector collimation 
thickness DCT is smaller than that of the single row 
radiation detector 13a shown in Fig. 5A, but images 
over a wider range can be taken at a time as a whole. 
The spatial resolution (body axis resolution) in the 
go-around axis direction of the tomographic image 



obtained improves as the detector collimation thickness 
becomes smaller. 

Next, the processing of creating a three- 
dimensional tomographic image of the object 12 by the 
reconfiguration means 22 from the projection data 
detected by the radiation source detector 13 will be 
explained . 

Fig. 6 is an example of flow chart showing 
the processing according to the Feldkamp 
reconfiguration method. 

On the other hand, Fig. 9 is a flow chart 
showing a processing operation of the reconfiguration 
means 22 of the tomograph according to an embodiment of 
the present invention. This flow chart assumes 
processing carried out slice by slice. 

The reconfiguration means 22 in Fig. 2 is 
provided with operating data phase range calculation 
means for determining a projection data phase range 
capable of back projecting for each reconfigured voxel, 
cone angle correction means for multiplying each row of 
projection data by a coefficient which is dependent on 
the angle of inclination of radiation from the 
radiation source, one-dimensional rearrangement 
processing means for obtaining parallel beam projection 
data from fan beam projection data obtained from a fan- 
shaped fan beam viewed from the go-around axis 
direction generated from the radiation source, filter 
correction means for superimposing the reconfiguration 
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filter on the parallel beam projection data and 
creating filter-processed parallel beam projection data 
and parallel beam three-dimensional back projection 
means for carrying out three-dimensional back 
projection on the filter-processed parallel beam 
projection data to the back projection region 
corresponding to a region in concern based on the 
determined projection data range capable of back 
projection. 

Based on the above described configuration, 
in Fig. 9, the operating data phase range calculation 
means determines the data range used for each voxel in 
step S4 first. Next, in step S5, the cone angle 
correction means multiplies each row of the projection 
data by a coefficient which is dependent on the angle 
of inclination of radiation and in step S6, the one- 
dimensional rearrangement processing means associates 
the fan beam projection data obtained from a fan-shaped 
fan beam viewed in the go-around axis direction 
generated from the radiation source with the parallel 
beam projection data. Then, in step S7, the filter 
correction means superimposes the reconfiguration 
filter on the parallel beam projection data and 
generates filter-processed parallel beam projection 
data. Next, in step 8, the parallel beam three- 
dimensional back projection means performs three- 
dimensional back projection on the filter-processed 
parallel beam projection data to the back projection 
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region corresponding to the region in concern based on 
the determined projection data range capable of back 
pro j ection . 

Next, the respective steps shown in Fig. 9 
will be explained. 

First, in step S4, the operating data phase 
range calculation means determines each data range used 
for each of all voxels in the slice. 

In the geometry shown in Fig. 28A, 28B and 
Fig. 29, suppose the distance between the radiation 
source 11 and rotation center is SOD, the relative 
movement distance in the body axis direction of the 
radiation source 11 with respect to the examinee per 
rotation of the scanner on the radiation detector 13 
(e.g., amount of feeding of the table) is J, the go- 
around phase of the fan beam source is P, beam 
spreading angle between the beam directed to the 
reconfigured voxel I (x, y, z) and central beam is a 
and the go-around phase of the parallel beam is <|>. 
Then, the fan beam source position S(p)=S(x s , y s , Zs) is 
expressed by the following Expression 2. Furthermore, 
when this is rearranged and replaced by a parallel 
beam, the fan beam source position is expressed by 
Expression 3. 
[Expression 2] 



S (j3) =S (SOD - sin^, -SOD-cosj3, J j3 / 2 tt ) 



[Expression 3] 
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S ( <t> ) =S (SOD- sin ( <j> + a ) , -SOD- cos (<*> + a),J(0 + a)/2 7r) 

Here, suppose the traveling direction of the 
parallel beam is W, the direction perpendicular to this 
traveling direction W (channel direction of parallel 
beam) is T. Then, the T coordinate and W coordinate 
5 when the parallel beam at phase (j> passes through 

coordinate (x, y) are expressed by Expression 4 and 
Expression 5 respectively. 
[Expression 4] 

T (x, y, (£>) =x-cos(£> + y-sin<£ 
[Expression 5] 

W(x, y, <£>)= — x-sin0-f-y-cos<f> 

Furthermore, the distance s_tz_dist between 
the X-ray source and T-Z plane (plane passing through 
10 the go-around axis and perpendicular to the parallel 
beam) is expressed by the following Expression 6. 
[Expression 6] 

s_t z_d i s t (x, y, <f>) = (SOD 2 -T (x, y, 0) 2 ) 1/2 
Furthermore, when the parallel beam with 
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phase (() passes through the reconfigured voxel I (x, y, 
z) and crosses the radiation detector 13 whose distance 
from the radiation source 11 is SID, suppose the 
coordinates of a system formed of the V axis (the same 
go-around axis direction as the z axis, the origin 
position thereof is detector center) of the radiation 
detector 13 and the X-Y axis are H (x, y, <J>) . Then, the 
coordinates are expressed by Expression 7. Here, while 
the Z axis matches the V axis, they are different in 
that the Z axis uses the scan start position as the 
origin position and the V axis uses the detector center 
as the origin position. 
[Expression 7] 

H (x, y, 0 ) = ( z- J ( <t> + a ) /2tz ) *SID/ (s_tz_dist (x, y, <f> ) +W (x, y , <f> ) ) 

Here, in Fig. 28A, oc=arcsin ( t /SOD) =A' t+B' , A' 
and B' are coefficients of approximate straight lines 
obtained by approximating arcsin ( t /SOD) . 

Based on Fig. 28C, how the phase range 
calculated in step S4 is determined will be explained.' 
As already described, this phase range differs from one 
voxel to another, and therefore the phase range is 
determined for each voxel and expressed by f here. The 
phase range is determined by selecting one which 
narrows the plane spreading angle of the cone beam with 
respect to the vertical line in the detector direction, 
that is, cone angle most. f is normally between 1 and 



2. At B s and B e which are the ends of phase range fn 
used to back project the reconfigured voxel I (x, y, 
z) , in order to minimize the absolute value of the cone 
angle, it is possible to select a go-around phase <J) of 
the parallel beam so that the difference in the 
absolute values between the coordinates H(x, y, §+fn) 
(corresponds to finally determined cone angle) and 
coordinates H(x, y, <J>) (corresponds to initial cone 
angle) become the smallest possible value. More 
specifically, an example of calculation algorithm of 
the data range with the smallest cone angle is shown in 
Fig. 10. As shown in step S101, when it is assumed 
that the initial value of <j) is -frc/2, calculation phase 
accuracy is Q (phase angle per view is normally used as 
Q, but when priority is given to the processing time, a 
phase angle exceeding one view can also be used) , the 
sum of the H(x, y, $+fn) and H(x, y, <|>) is err(x, y, <|>) 
(hereinafter expressed as "err") and a minimum value of 
this sum err is err_min (initial value is err_min=f n) , 
err is expressed by Expression 8 and Expression 9 shown 
in step, S102 . 
[Expression 8] 

e r r =H (x, y, <f> ) +H (x, y, 0 + f tz ) 
[Expression 9] 

if [err m i n > e r r ] , err m i n = e r r 
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Here, when <J> increases , err decreases and when 
<J> decreases, err increases, and so the following 
Expression 10 and Expression 11 are repeated. 
[Expression 10] 

if [err>0], 0 = 0 +Q 

[Expression 11] 

if [ e r r < 0] , <t> = $ —Q 

Through this repeating processing, if err is 
5 compared with err_min as shown in step S103, when err 
becomes a minimum value, minimum values appear 
repeatedly and err=err_min, and so by carrying out 
repeating processing until err=err_min is obtained, it 
is possible to select § as shown in step S104 so that 

10 the difference in the absolute values between 

coordinates H(x, y, <j>+f7u) and coordinates H(x, y, (J>) 
becomes the smallest possible value as shown in step 
S105. If the decision in step S105 results in err>0, 
<j>=<t>+Q as in step S107, and if err<0, <(>=<|)-Q as in step 

15 S108. Thus, the phase range (Bs<(j><Be) is expressed by 
the following Expression 12 and Expression 13. 
[Expression 12] 

B s ( x , y , z ) = 0 



[Expression 13] 
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Be (x, y, z) = <f> + f 7c 

Here, the phase range has been determined 
using the simplest method as described above, but this 
is the problem of calculation of a minimum value of the 
function err ((f)) in the phase range (-£n/2<$<fn/2) and 
it is also possible to use an existing method, for 
example, Brent's method and golden division method 
(golden section search) and combine various methods to 
calculate <j> and ty+fn so that err((j>) becomes a minimum 
value. Furthermore, it is also possible to increase 
the processing speed by optimizing the initial value 
when <(>=-f7r/2 is determined. 

Furthermore, with regard to the phase range 
fTC used to back project the reconfigured voxel (x, y, 
z) , it is also possible to determine Bs and Be by 
determining go-around phase <|> of the parallel beam so 
that the absolute value of the angle of inclination of 
the beam (cone angle) of the X-ray beam becomes small 
at the end of the phase ranges and extend the data 
range to both ends of the data range as shown the 
following Expression 14 and Expression 15. 
[Expression 14] 



Bs (x, y, z) =d>~ (f-1) 7i/2 



[Expression 15] 



Be (x, y, z) = <t> + f tu + (f-1) n / 2 

Next, the cone angle correction step of 
multiplying each row of the parallel beam projection 
data by a coefficient which is dependent on the cone 
angle using the cone angle correction means in step S5 
shown in Fig. 9 will be explained. 

Filter correction in reconfiguration is 
filtering corresponding to the distance from the go- 
around axis in the reconfigured image and it is 
necessary to apply a filter corresponding to the cone 
angle to correct the influences of beam inclination. 
Here, suppose data before filter correction is P par a(<f>, 
t, v) , data after filter correction is fP para {§, t, v) , 
and the reconfiguration filter function is g(t). Then, 
the reconfiguration filter processing can be expressed 
as shown in Expression 16 using a convolution method 
and of this, the cone angle correction is the portion 
expressed by Expression 17. As is evident from 
Expression 16, since the cone angle correction term is 
a coefficient corresponding to the detector row 
position v (cone angle) , cone angle correction can be 
carried out both before and after filter correction. 
For this cone angle correction, a publicly known 
technology in the three-dimensional back projection 
techniques including a Feldkamp method is applicable. 
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[Expression 16] 

fPparaU> t, v) = | S ID P p . r a ( <fr , t - t \ v) g ( t ') d t ' 

J ^°VS I D 2 + v 2 

where t' is a variable of integration in Expression 16. 
[Expression 17] 

S I D/Js I D 2 +v 2 

Next, the rearrangement processing 
(rebinning) using the one-dimensional rearrangement 
processing means in step S6 shown in Fig. 9 will be 
explained. Here, Fig. 32A and Fig. 32B show a 
relationship between a fan beam and a parallel beam. 
Figs. 32A to C show a 180° reconfiguration of a fan beam 
and Fig. 32D shows a 180° reconfiguration of a parallel 
beam. When X-ray beams (SI to S3) irradiated in the 
same vector direction viewed from the go-around axis 
direction are gathered together, it is possible to 
virtually create a parallel beam as shown in Fig. 32B. 

In order to enhance the calculation speed, 
one-dimensional rearrangement processing is carried out 
which rearranges a fan beam irradiated in a fan shape 
viewed from the go-around axis direction as shown in 
Fig. 11A and Fig. 11B into parallel beams which are 
parallel viewed from the go-around axis as shown in 
Fig. 12A (Fig. 12B is an exploded view of Fig. 12A) and 
Fig. 12B. Furthermore, Fig. 13A and 13B show the 
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parallel beams rearranged in the go-around axis 
direction, which will be described later. 

Figs. 11B, 12B and 13B are exploded views of 
beams and their respective focuses on the detector 
corresponding to Figs. 11A, 12A and 13A. If the fan 
beam is represented by P fan (P, a, v) and parallel beam 
is represented by P par a(<t>, t, v) , the fan beam spreading 
angle in the go-around direction ct=arcsin ( t/SOD) and 
p=(|>+a (see Fig. 28A) , and therefore the rearrangement 
processing can be expressed by Expression 18. 
[Expression 18] 

P para ( <t> > t, v) =P fan ( + a , a , v ) 

Next, a convolutional calculation (filter 
correction processing) of a reconfiguration filter 
carried out to correct blurs of projection data using 
the filter correction means in step S7 shown in Fig. 9 
will be explained. 

For filter correction, two types of methods; 
a convolution method which carries out a convolutional 
calculation in a real space and a Fourier method which 
carries out a multiplication in a Fourier space. The 
convolution method in the former is convolutional 
processing on a filter function which has been inverse 
Fourier transformed in a real space. The Fourier 
method in the latter is processing consisting of 
transforming into a Fourier space using a Fourier 



transform, multiplying it by a filter function (spatial 
frequency filter) and then applying an inverse Fourier 
transform. 

Both are processes mathematically equivalent, 
but filter processing in a Fourier space which provides 
high-speed calculation is generally used. For the 
filter used for reconfiguration, it is possible to 
select and use Shepp and Logan, Ramachandran and 
Lakshminarayanan or these filter functions which have 
been modified through clinical experiences based on 
clinical experiences. Suppose the parallel projection 
data is P pa ra(<|>/ t, v) , the parallel projection data 
after filter processing is fP par a(<K t, v) and the 
reconfiguration filter is G ( co ) . Then, Fourier space 
filtering according to a Fourier method can be 
expressed by Expression 19. 
[Expression 19] 

fp P ara(^ > v ) = ~ I ( <t> , t, v) • exp ( — i co t) dt • G ( co ) • exp (i co t) dco 

4 7C J -00 J-00 

On the other hand, when an inverse Fourier 
transform g(t) of G ( co ) is expressed as shown in 
Expression 28, the real space filtering according to 
the convolution method can be expressed by Expression 
29. 

[Expression 20 ] 
1 f 00 

g ( t ) = - — G(oo) -exp (icot) dco 

Z 71 J-00 
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[Expression 21] 

~ oo 

fp psr ,(i t, v) = I fo, t-t\ v) g ( t ') d t' 

J— OO 

where t' is a variable of integration in Expression 21. 

For simplicity, the direction in which the 
filter is applied is assumed to be the T direction 
here, but it is possible to apply the filter in a high- 
dimensional direction combining the V direction, T 
direction and (J) direction. Furthermore, the projection 
data is handled as continuous data here, but since the 
projection data is actually discrete data, it is 
necessary to use a publicly known interpolation method 
to calculate the projection data in a discrete manner. 
This discrete calculation method has been practiced so 
far and is similar to filter correction, etc., used for 
weighting spiral correction reconfiguration. 

Further, implementation of the three- 
dimensional back projection corresponding to the data 
range determined by the aforementioned determining 
means in step S8 of Fig. 9 will be explained. 

As shown in Figs. 28A and 28B, if the 
reconfigured voxel is I(x, y, z) , the V axis direction 
position that matches the go-around axis on a 
cylindrical detector centered on the radiation source 
11 is v and the position on the T axis substantially 
orthogonal to this V axis is t, then the reconfigured 
voxel I(x, y, z) is expressed by Expression 22. 
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[Expression 23] 



I ( x , y , z) =- [ 



Be (x, y , z ) 



7C J 



B s (x, y , z) 



f P 



para 



((f) t t , V ) d 0 



where t = xcos<|> + ysin<j> 




SODcosa — xsin<f> + ycos<£ 2 71 



SID , T 



This algorithm handles the projection data 
and reconfigured image which should originally be 
handled discretely as continuous data, and therefore it 
is actually desirable to use an interpolation method 
such as Lagrange interpolation and calculate discretely 
through interpolation in three directions of the phase 
direction, detector row direction and detector channel 
direction. To realize a high-speed calculation at the 
sacrifice of accuracy, the above described v can also 
be v= ( z-J(j>/27i) • SID/ (SODcosa - xsin(f> + ycos<|>) . 



dimensional back projection method, it is possible to 
obtain an image of good quality with fewer errors 
compared to the conventional two-dimensional 
reconfiguration (weighting spiral correction method) . 
Furthermore, by performing back projection from data 
with minimum errors (data with a small cone angle) for 
each voxel to achieve higher image quality and 
including determining means for determining the data 



According to such a filter correction three- 
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phase range used for each voxel to realize this, or 
more specifically, determining the data phase range for 
each voxel so that the absolute values of the angles of 
inclination of radiation beams become the same at both 
ends of the data in the phase range, it is possible to 
use projection data with a smaller cone angle and by 
performing corrections using a weighting function for 
each voxel while maintaining redundancy, it is possible 
to obtain an image with discontinuity in the data phase 
direction reduced. 

Especially, by using data 270 degrees in the 
phase direction and performing weighting as shown in 
Fig. 20A, it is possible to reduce discontinuity at the 
ends of data to a minimum. With this 270-degree data, 
it is possible to correct the discontinuity at 180- 
degree data ends as shown in Figs. 14A and 14B using 
the data phase with the least discontinuity with a 
phase difference of 90 degrees as shown in Figs. 15A 
and 15B. Data becomes discontinuous at the position n 
in Figs. 14A and 14B, which may cause artifact. That 
is, it is possible to reduce the data discontinuity to 
a minimum and realize reconfiguration of higher 
quality. Further, if it is possible to accurately 
equalize absolute values of the angle of inclination of 
the radiation beam (cone angle) at both ends of data, 
it is also possible to calculate the detector row 
direction position from the data start direction and 
end direction simultaneously and realize a high-speed 



calculation. Furthermore, since the same phase range 
is used during back projection of each reconfigured 
voxel, a weighting function for redundancy correction 
is determined by a single expression, which allows a 
high-speed calculation. Errors are mutually corrected 
at positions of ti/2, n in Fig. 15B. 
Embodiment 2 

Fig. 16 is a flow chart showing a processing 
operation of the reconfiguration means 22 according to 
another embodiment of the present invention. 

In this embodiment, a filter correction is 
performed in step S7 shown in Fig. 9, and then 
redundancy correction weighting is performed in step S9 
and three-dimensional back projection is performed in 
step S8. 

In the case of this embodiment, the 
reconfiguration means further includes redundancy 
correction weighting means for realizing a redundancy 
correction using a weighting function whose shape 
changes according to the phase width with respect to 
filter-processed projection data over a projection data 
range fn determined and obtained by the operating data 
phase range means. 

Processes in steps S4 to S8 are the same as 
the procedure already explained using Fig. 9 and 
therefore a weighting process using the redundancy 
correction weighting means in step S9 will be explained 
here . 
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As shown in Figs. 17 to 19, data of 180 
degrees or more is used for each voxel to reconfigure 
an image and data correction is performed through 
weighting using the weighting function as shown in Fig. 
5 20A to correct data redundancy. More specifically, as 
shown in the weighting function W(0) in Figs. 21A and 
21B and Expression 23 to Expression 25, weighting is 
performed on the phase data range which varies from one 
voxel to another so that the sum of weights at the same 

10 phase and opposite phase used for back projection 

remains equal at the respective phases. Here, when the 
data width used for each voxel is B=fn, when B=tt (when 
f=l) , weighting becomes as shown in Fig. 21A, likewise 
when B=3n/2 (when f=3/2), weighting becomes as shown in 

15 Fig. 21B and when B=2tt (when f=2), weighting becomes as 
shown in Fig. 21C. 
[Expression 23 ] 

W (0) = ( (B/2) +6) /B-K 
where [-n/2<9< (2tc-B) /2] . 
[Expression 24] 
W (0) = 1 



where [- (2tt-B) /2<9<(27i-B) /2] . 
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[Expression 25] 

W (0) = ( (B/2) -6) /B-7u 

where [ (2tc-B) /2<9<B/2] . 

In Fig. 20A and Figs. 21A to 21C, weighting 
which linearly changes in the view direction is 
performed, but it is also possible to perform weighting 
which changes nonlinearly in the view direction as 
5 shown in Fig. 20B. The nonlinear weighting function 
W (0) shown in Fig. 20B can be calculated, for example, 
from the above described weighting function W(6) as 
shown in Expression 26 to Expression 28. Furthermore, 
only the case with B < 2n is described here, but the 
10 case with B>27i can also be easily calculated based on a 
similar concept . 
[Expression 26] 

W (9) =3 (W (0) ) 2 -2 (W (9) ) 3 
where [ -n/2<Q< ( 2tt-B) /2] . 
[Expression 27] 
W (0) - 1 



where [- (2tt-B) /2<9</ (2tc-B) /2] . 



[Expression 28] 



W (0) =-3 (W (0) ) 2 +2 (W (0) ) 3 

where [ (2rc-B) /2<9<B/2] . 

In the above described phase range 
calculation process for each voxel, such a tomograph is 
a three-dimensional reconfiguration method which 
determines a phase range of fit [rad] in the view 
direction and carries out a redundancy correction using 
a weighting function, and provides data with redundancy 
(extends the back projection phase width beyond 180 
degrees) , assigns weights using the weighting function, 
and can thereby reduce discontinuity at the data ends 
(at the start/end of image taking) and obtain an image 
with the influence of movement of the examinee reduced 
to a minimum. 

When a fan beam reconfiguration is used in a 
Wang method or IHCB method of the conventional 
examples, the redundancy (projection phase range) of 
projection data obtained varies from one voxel to 
another. For example, when the radiation source 
performs back projection from data obtained by rotating 
the phase by 180 degrees as shown in Fig. 22, the data 
phase range capable of back projection varies from one 
reconfiguration pixel to another and data in a phase 
range of 180 degrees or more is obtained at pixel a, 
but only data of 180 degrees or less is obtained at 



pixel b. Thus, because data redundancy varies from one 
pixel to another, when back projection is performed 
from projection data of 360 degrees or less, 
complicated redundancy correction processing is 
required at the time of back projection. In a three- 
dimensional reconfiguration in particular, a cone angle 
needs to be considered, and therefore more complicated 
redundancy correction processing is required, which 
constitutes one of causes of an increase in the 
calculation time. Furthermore, this redundancy 
correction processing is also associated with measuring 
throughput (relative moving speed between the focus and 
examinee) . Unlike these conventional examples, the 
present application rather takes advantage of 
redundancy, uses data with a phase range of 180 degrees 
or more for each voxel and thereby prevents generation 
of discontinuity due to movement, etc., and also 
improves the data efficiency. 
Embodiment 3 

Fig. 23 is a flow chart showing a processing 
operation of reconfiguration means 22 according to a 
further embodiment of the present invention. 

As shown in Fig. 23, after processes in step 
S4 and step S5, this embodiment carries out a 
rearrangement process on projection data of an image 
taken with a view of a multiple of 4 in step Sll. 
Then, in step S7, a filter correction is performed and 
in step S12, projection data whose phase in the go- 



around direction differs by Ntt/2 (N=l, 2, 3, ) 

[rad] is grouped by grouping means and in step S9, a 
redundancy correction weighting process is carried out 
and in step S8, the grouped projection data is back 
projected to a square image group by group. 

In order to realize such processing, means 
for acquiring projection data whose number of images 
taken per rotation is a multiple of 4 is provided, the 
reconfiguration means 22 includes means for 
superimposing a filter on this projection data, 
grouping means for grouping data at the same channel 
position and whose projection phase in the go-around 

direction differs by Ntt/2 (N=l, 2, 3, ) [rad] 

and back projection means for back projecting into a 
square image array group by group using this grouping 
means . 

Thus, in order to enhance the speed of back 
projection which takes a maximum calculation time in 
creating an image, an image is taken with a view of a 
multiple of 4 and a fan beam is reconfigured in Fig. 

23, and data is converted to data whose number of views 
is a multiple of 4 through the rearrangement process in 
step S13 and a parallel beam is reconfigured in Fig. 

24, taking advantage of the fact that the reconfigured 
image array is square and images are taken while the 
detector is going around the reconfigured image. 

In both cases, projection data whose phase in 
the go-around direction differs by Ntt/2 (N=l, 2, 



3/ ) [rad] is grouped and back projection is 

performed on the square image group by group, and 
therefore it is possible to reduce, for example, the 
number of calculations of the channel direction 
position in a full reconfiguration and interpolation 
coefficient to 1/4 (it is possible to reduce the number 
of calculations to 1/2 in a half reconfiguration) . 
This is because if the reconfigured image is square, 

data whose phase differs by Ntt/2 (N=l, 2, 3, ) 

[rad] and the square which is the reconfigured image 
have the same positional' relationship. 

Furthermore, the number of views is set to a 
multiple of 4 is to accurately calculate data whose 

phase differs by Nti/2 (N=l, 2, 3, ) [rad]. 

Furthermore, in both cases of full reconfiguration and 
half reconfiguration, it is possible to create images 
by calculating the channel position within a range of 
1/4 (7i/2[rad]) of one revolution. In terms of a full 
reconfiguration, the amount of calculation becomes 1/4 
and though the calculation is performed using only one 
calculator, it is possible to obtain a result close to 
the case where parallel calculations are performed 
using four calculators.. That is, it is possible to 
realize high performance at a low cost. Needless to 
say, it is also possible to set the number of views to 
a multiple of 4 during image taking and perform 
reconfiguration directly from a fan beam without any 
rearrangement process (rebinning) . Furthermore, when a 
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display pixel is a hexagon, it is possible to group 
projection data whose phase in the go-around direction 

differs by Ntt/3 [rad] (N=l, 2, 3, ) and back 

project to the hexagonal image group by group. When 
the display pixel is polygonal and has C sides, the 
above described phase in the go-around direction is 
2tt/C [rad] . 

Next, group-by-group back projection will be 

explained. 

First, as shown in Fig. 25, when only the X-Y 
plane is considered and a beam passing through a 
voxel (x, y) irradiated from a focus position S(p) 
having phase P is irradiated to a position u on the 
detector, group-by-group back projection processing is 
expressed by Expression 29 to Expression 32. 



calculation in the channel direction is 1/4 (view=4N; N 
is an integer) and the calculation start position is 
not limited. Furthermore, reference numeral 131 
denotes a reconfiguration region. 
[Expression 29] 



Fig. 25 shows a case where the amount of 
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f P 



para 
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[Expression 30] 
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[Expression 31] 
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(-x, -y) =i f 

71 J 



B s (x, y) 
B s (x, y) 



71 

2 



fP para(* + ^' t > v ) d <t> 



[Expression 32] 

I (-y, X) = - I , ■ 2 f Pp a r a t, v d 0 

7Z JB s (x, y) L 2 J 

A beam irradiated from phase P+ti/2 and passing 
through a voxel (-y, x) is irradiated to the position u 
on the radiation detector as in the case of being 
irradiated from phase p to a voxel (x, y) . Likewise, a 
5 beam irradiated from phase p+7t passes through a voxel (- 
x, -y) and is irradiated to the position u on the 
radiation detector. Likewise, a beam irradiated from 
phase P+3tt/2 passes through a voxel (y, -x) and 
irradiated to the position u on the radiation detector. 

10 Thus, by performing back projection from the grouped 
data to four pixels which use the same radiation 
detector position data, it is possible to calculate the 
radiation detector position and reduce the number of 
times interpolation parameters are calculated. 

15 Fig. 26 shows a case where the amount of 

calculation in the row direction is l/2n (where n 
denotes the number of data revolutions) and rpitch=J/2N 
(where N is an integer) . As shown here, suppose an x- 
y-z space (Euclidean space) is considered and a 
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relative moving speed between an object and radiation 
source in the go-around axis direction (e.g., bed 
feeding speed) is J and a beam passing through a voxel 
I (x, y, z) irradiated from a focus position S(P) which 
is phase p is irradiated to the position v in the go- 
around axis direction on the radiation detector. A 
beam irradiated from phase P+271 and passing through a 
voxel (x, y, z + J) is irradiated to the position v in the 
go-around axis direction on the radiation detector as 
in the case where the beam is irradiated from phase P 
to the voxel I (x, y, z) . Likewise, a beam irradiated 
from phase P+tt passes through a voxel I (-X, -y, z+J/2) 
and is irradiated to the position v in the go-around 
axis direction on the radiation detector. Talcing 
advantage of this, the object is associated with the 
relative moving speed of the radiation source in the 
go-around axis direction according to the 
reconfiguration intervals and data pieces with phases 

differing by Ntt (N=l, 2, 3 ) [rad] from one 

another are grouped and back projected group by group. 

According to such group-by-group back 
projection, by associating the pixel intervals of the 
voxel in the body axis direction at MDCT with the 
object and the relative moving speed of the radiation 
source in the go-around axis direction, it is possible 
to calculate the position in the body axis direction at 
high speed, and when an image is created from data of a 
plurality of revolutions obtained by taking images 



through a spiral scan, it is possible to enhance the 
speed of back projection which takes a maximum time to 
create images. 

Here, the spiral period in the body axis 
direction is synchronized with the period of the 
reconfigured voxel in the body axis direction, and 
when, for example, the pixel interval (voxel pitch) in 
the body axis direction is rpitch [mm] , the relative 
moving speed (bed moving speed) of the radiation source 
in the body axis direction with respect to the examinee 

is set to table=2-N- rpitch (N=l, 2, 3, ). In 

this way, at the phase of the radiation source which is 

N71 (N=l, 2, 3, ) [rad] , the positions on the 

radiation detector at which the beams passing through 
the voxel I (x, y, z) whose body axis direction 
position is Z [mm] and the voxel I (-x, -y, N-J/2+Z) 
whose body axis direction position is (N-J/2)+Z [mm] 
intersect with each other are the same, and therefore 
calculating a beam passing through a voxel with a view 
at the time of back projection is equivalent to 
simultaneously calculating the row positions at phases 

differing by Nn (N=l, 2, 3, ) [rad] from each 

other. Thus, calculations of the row direction 
positions of the radiation detector and interpolation 
coefficients over the total measuring range are 
completed within the7i [rad] range in the view 
direction . 

In the above described embodiment, no 
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rearrangement in the radiation detector row direction 
is performed so that descriptions in the rearrangement 
processing do not become complicated, but to enhance 
the speed of back projection, it is also possible to 
perform rearrangement in the row direction on the plane 
located at the rotation center which crosses the 
parallel beam at right angles as expressed in P par a(P, t, 
v) =P fa n (<f>+a, a, (SID/SOD-cos (a) )• (v-J-a/27t) ) where 
a=arcsin (t/SOD) as shown in Fig. 13A so that points of 
intersection in the parallel beam channel direction 
become the same v coordinates. When such a 
rearrangement in the row direction is performed, it is 
possible to reduce the number of arcsin calculations 
used to calculate a at the time of back projection and 
realize faster processing. In this case, the operating 
data phase range for each voxel can be likewise 
calculated with H(x, y, <|>) in the above described 
expression changed to Expression 33. 
[Expression 33] 

H (x, y, <*>) = fz-*—± s_t z_d 1st (x, y, <t>) 

^ 2tt J s tz dist ( x , y, </> ) + w ( x , y x ) 

Furthermore, in this case, v in Expression 22 
is changed to v= ( z- ( J/2n) (<|>+a) ) • SODcosa/ (SODcosa - 
xsin<|>+. ycos(|)) to obtain the projection beam used for 
back projection. 

Furthermore, the tomograph in the above 
described embodiment is also applicable to products 
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using X-rays, gamma rays, neutron rays, positron, 
electromagnetic energy or light. Furthermore, the scan 
system is not limited to any of first-generation to 
fourth-generation systems and this tomograph can also 
be used for a multi-tube CT incorporating a plurality 
of radiation sources and doughnut type tube CT . 
Furthermore, with regard to the shape of the radiation 
detector, this tomograph is also applicable to any 
radiation detector such as detectors arranged on a 
cylindrical surface centered on the radiation source, 
plane detectors,, detectors arranged on a spherical 
surface centered on the radiation source and detectors 
arranged on a cylindrical surface centered on a go- 
around axis, etc. Furthermore, the position of a 
radiation detector corresponding to the reconfigured 
voxel is calculated every time, but for grouping in the 
channel direction, it is also possible to store a table 
of reconfiguration parameters calculated beforehand 
corresponding to N/4 revolutions (0<P<N7t/2, N=l, 2, 

3 ) in a memory, read this stored parameter table 

at the time of reconfiguration and realize 
reconfiguration based on this parameter table. 
Adopting such a configuration allows calculations of 
addresses corresponding to 4 views all at once. Such 
simplification of calculations is a technique 
unparalleled in conventional examples. The above 
described N/4 revolutions apply to the case where the 
shape of a display pixel is rectangular and when the 



display pixel is hexagonal, data can also be grouped 
every N/6 revolutions. 
Embodiment 4 

Fig. 27 is a flow chart showing a processing 
operation of reconfiguration means 22 in a tomograph 
according to an embodiment of the present invention. 

First, the reconfiguration means 22 is 
provided with operating data phase range calculation 
means for determining a projection data phase range 
capable of back projection for each reconfigured voxel, 
approximate straight line calculation means for 
calculating an approximate straight line for a curve 
indicating the radiation source position with respect 
to the channel direction position corresponding to a 
region in concern of parallel beam projection data 
obtained by a parallel beam of a parallel shape viewed 
from the go-around axis direction generated from the 
radiation source, cone angle correction means for 
multiplying each row of projection data by a 
coefficient which is dependent on the angle of 
inclination of radiation from the radiation source, 
one-dimensional rearrangement processing means for 
obtaining parallel beam projection data from the fan 
beam projection data obtained by a fan-shaped fan beam 
viewed from the go-around axis direction generated from 
the radiation source, filter correction means for 
superimposing a reconfiguration filter on the parallel 
beam projection data and creating filter-processed 
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parallel beam projection data and parallel beam three- 
dimensional back projection means for three-dimension 
back projecting the filter-processed parallel beam 
projection data to a back projection region 
corresponding to a region in concern along the 
approximate irradiation trace of the radiation beam 
calculated using the approximate straight line based on 
the determined projection data range capable of back 
projection. 

Based on the above described structure, the 
data range used for each voxel is determined using the 
operating data phase range calculation means in step S4 
first, and an approximate straight line for a curve 
indicating the radiation source position with respect 
to the channel direction position of the parallel beam 
projection data obtained by a parallel beam of a 
parallel shape viewed from the go-around axis direction 
generated from the radiation source by the approximate 
straight line calculation means is calculated in step 
S14. Next, in step S5, the cone angle correction means 
multiplies each row of the projection data by a 
coefficient which is dependent on the angle of 
inclination of radiation and in step S6, the one- 
dimensional rearrangement processing means associates 
the fan beam projection data obtained from a fan-shaped 
fan beam viewed from the go-around axis direction 
generated from the radiation source with the parallel 
beam projection data. Then, in step S7, the filter 
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correction means superimposes a reconfiguration filter 
on the parallel beam projection data and creates 
parallel beam projection data subjected to filter 
processing- Then, in step 15, based on the projection 
data range capable of back projection determined by the 
parallel beam three-dimensional back projection means, 
the parallel beam projection data subjected to filter 
processing is three-dimensional back projected to the 
back projection region corresponding to the region in 
concern along the approximate irradiation trace of the 
radiation beam calculated using an approximate straight 
line . 

Steps S4 to S7 are the same as those already 
explained in other embodiments. 

The calculation of an approximate straight 
line by the approximate straight line calculation means 
in step S14 for the curve indicating the radiation 
source position with respect to the channel direction 
position of the parallel beam projection data obtained 
by a parallel beam of a parallel shape viewed from the 
go-around axis direction generated from the radiation 
source will be explained. 

Here, a technique using a least squares 
method will be shown. First, when an approximated 
curve and an approximate curve will be considered. A 
coordinate Z ± of the focus at the channel i position of 
a parallel beam is expressed by the following 
Expression 34 and an approximate straight line z A except 
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an arcsin calculation is expressed by the following 
Expression 35. Here, suppose the position of the 
channel i of the parallel beam in the t axis direction 
is ti . 

[Expression 34] 

z i = J • arcs in (tj/SOD) / 2 n 

[Expression 35] 

z A ( t i ) =A • t i+B 

A, B in the expression can be calculated more 
specifically as follows. 

When points on the approximated curve within 
a diameter FOV of the circular region in concern shown 
in Figs. 28A and 28B are approximated with an 
approximate straight line and an evaluation function 
for minimizing the error between the approximated curve 
and approximate straight line is used according to a 
least squares method, the points can be expressed by 
Expression '36. N t denotes the number of samples. 
[Expression 36] 

E 2 (A, B) = ]T( Zi - z A ( ) 2 =|](z i -A- ti-B) 2 
i = 1 i = 1 

Here, Expression 36 is minimized to determine 
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A, B. With the minimum values, the differentiation 
values with respect to A and B of Expression 44 are 
zeros as shown in Expression 37 and Expression 38. 



[Expression 37] 

0=^=-2£(z-A. t-B) 
a& i = i 

[Expression 38] 

*F 2 N i 
0 = ^=-2l(V (z-A- t-B) } 

i = 1 

For simplicity, when the following sums in 
5 Expression 39 are introduced and these sums are 

substituted into Expression 36 and Expression 37, then 
Expression 40 and Expression 41 are obtained. 
[Expression 39] 




[Expression 40] 



B • S +A • S t = S z 



[Expression 41] 



B • S t +A • S tn = S t z 



The solution of these simultaneous equations 
is given in the following Expression 42 to Expression 
44 . 

[Expression 42] 

A = S • S t t - (S t ) 2 

[Expression 43] 

A= (S t t • S z -S t • S t z ) /A 

[Expression 44] 

B= (S • S tz -S t • S z ) /A 

Thus, by ■ substituting this into z A ( t± ) =A • t±+B 
shown in Expression 35, it is possible to obtain 
Expression 45 . 
[Expression 45] 

z A ( t j) = ((S tt • S 2 -S t • S t ,)/A)-t, + (S • S tz -S t - S 2 )/A 

Next, based on the determined projection data 
range capable of back projection in step S15 shown in 
Fig. 27, the parallel beam three-dimensional back 
projection means which performs three-dimensional back 
projection on the parallel beam projection data 
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subjected to filter processing to the back projection 
region corresponding to the region in concern along the 
approximate irradiation trace of the radiation beam 
calculated using an approximate straight line will be 
explained . 



reconfigured voxel is I (x, y, z), the relative 
movement distance of the radiation source 11 with 
respect to the examinee per rotation of the scanner on 
the radiation detector is J, the go-around axis 
direction position on the cylindrical radiation 
detector 13 centered on the radiation source 11 is v, 
the position on the T axis substantially perpendicular 
thereto is t and the coordinates are T(x, y, , then 
Expression 46 to Expression 50 are obtained 
respectively. In Fig. 29, reference character A 
denotes an image array of I (x, y, z) and 111 denotes an 
X-ray beam. 



dimensional back projection along the beam trace over 
the back projection data range determined by the data 
phase range calculation means. 



As shown in Fig. 28A and Fig. 29, suppose the 



Expression 46 shows a weighting three- 



Expression 50 shows a radiation beam trace 



calculated using an approximate straight line. 



[Expression 46] 



B e (x,y,z) 



r 




B s (x,y,z) 



fPparaU>t,v) • W 0 -B s (x f y,z) - 



f 7C 



2 J 



• d<t> 
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[Expression 47] 

L ( x , y , 4> ) =Vs OD 2 — t 2 — x*sin<2> + y-cos(f> 

[Expression 48] 

t (x, y, <j))=x-cos<J>+y-sin<|> 
[Expression 49] 

v= ( Zl -z s ) 'SID/L x, y) 

[Expression 50] 

J* 0 + arcsin — - — T 

Here, in the three-dimensional back 
projection, projection data and a reconfigured image 
which should actually be handled discretely are handled 
as continuous data, and therefore it is actually 
5 necessary to calculate the data discretely using a 

combination of interpolation in three directions of the 
phase direction (time direction), radiation detector 
row direction and radiation detector channel direction 
using a publicly known interpolation method such as 
10 Lagrange interpolation. 

As is evident from the above described 
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reconfiguration method, Expression 50 has a large 
calculation load with an arcsin calculation included in 
the calculation of the focus z position of the 
conventional parallel beam as seen from a comparison 
with Expression 1, but this arcsin calculation is 
replaced by an approximate straight line, and therefore 
it is possible to simplify an amount of calculation of 
the parallel beam three-dimensional back projection 
method and drastically enhance the processing speed. 

However, this reconfiguration method may 
involve the risk of deterioration of accuracy due to 
the use of the approximate straight line, but this 
error remains at such a level that even when the 
diameter of FOV of a circular region in concern is 410 
[mm] , the distance SOD between the focus and go-around 
axis is 600 [mm], the distance SID between the focus 
and detector is 1000 [mm] , the number of detector rows 
row is 64 [rows], the detector element direction size 
dapp is 1 [mm] and the relative moving speed T is 60 
[mm/rot], a maximum error is on the order of 0.023 [mm] 
and absolute error average is. on the order of 0.014 
[mm] . This error is an error on the order of 2% 
(maximum 4%) considering the measuring accuracy and the 
z direction width of the beam at the rotation center of 
0.6 [mm] and is at a totally insignificant level taking 
into consideration that noise is included in measuring 
data. That is, the approximate calculation will not 
lead to deterioration of image quality. 
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Furthermore, in the process of determining 
the phase range for each voxel shown in step S4, a 
phase range of fn [rad] is determined in the view 
direction and a three-dimensional reconfiguration 
method whereby redundancy correction is performed using 
a weighting function is used, and therefore by 
providing the data with redundancy (extending the back 
projection phase width beyond 180 degrees) and 
assigning weights using a weighting function, it is 
possible to reduce discontinuity at the data ends, that 
is, at the of start and end of image taking and obtain 
an image with the influence of movement of the examinee 
suppressed to a minimum. 

Furthermore, when a fan beam is rearranged to 
parallel beams and then one-slice reconfigured image is 
reconfigured through three-dimensional back projection, 
the conventional art uses the same back projection 
phase range for all voxels, and the z direction 
positions of the focuses of parallel beams are not the 
same in the channel direction, and therefore a maximum 
cone angle back projected at each voxel increases. 
That the maximum cone angle used increases means that a 
wider detector is required depending on the go-around 
axis z direction, that is, the relative moving speed in 
the z direction between the examinee and focus 
decreases and the measuring throughput deteriorates. 
However, in this embodiment, the maximum cone angle of 
the beam used for back projection is reduced as 
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described above, and therefore it is possible to 
reconfigure a detector which is narrow in the z 
direction and improve the measuring throughput. 
Embodiment 5 

Fig. 30 is a flow chart showing a processing 
operation of reconfiguration means 22 of a tomograph 
according to a still further embodiment of the present 
invention . 

Here, the reconfiguration means 22 consists 
of operating data phase range calculation means for 
determining projection data phase range capable of back 
projection for each reconfigured voxel, approximate 
straight line calculation means for calculating an 
approximate straight line for a curve indicating the 
radiation source position with respect to the channel 
direction position of parallel beam projection data 
obtained by a parallel beam of a parallel shape viewed 
from the go-around axis direction generated from a 
radiation source, cone angle correction means for 
multiplying each row of projection data by a 
coefficient which is dependent on the angle of 
inclination of radiation, one-dimensional rearrangement 
processing means for associating fan beam projection 
data obtained from a fan beam of a fan shape viewed 
from the go-around axis direction generated from the 
radiation source with the parallel beam projection 
data, filter correction means for superimposing a 
reconfiguration filter on the corrected projection data 
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and creating filter-processed projection data, 
redundancy correction weighting means for carrying out 
a redundancy correction on the filter-processed 
projection data over a projection data range fn 
determined by the operating data phase range 
calculation means using a weighting function whose 
shape changes according to the phase width and parallel 
beam three-dimensional back projection means for 
performing three-dimensional back projection to a back 
projection region along an approximate irradiation 
trace of the radiation beam calculated based on the 
approximate straight line obtained by the approximate 
straight line calculation means while carrying out 
weighting processing on the filter-processed projection 
data using this redundancy correction weighting means. 

As in the case of Fig. 27, such 
reconfiguration means 22 superimposes a reconfiguration 
filter on the parallel beam projection data using the 
filter correction means in step S7, generates filter- 
processed parallel beam projection data and then in 
step S9, carries out a redundancy correction on the 
filter-processed projection data created by the filter 
correction means using a weighting function by the 
redundancy correction weighting means over the data 
range fn determined by the operating data phase range 
calculation means. Then, while performing weighting 
processing using this redundancy correction weighting 
means, in step 15, the filter-processed parallel beam 



projection data is three-dimension back projected to a 
back projection region corresponding to a region in 
concern along an approximate irradiation trace of the 
radiation beam calculated using an approximate straight 
line based on the projection data range capable of back 
projection determined by the parallel beam three- 
dimensional back projection means. 

Details of each step have already been 
explained with the same step numbers assigned, and 
therefore explanations thereof will be omitted. 

The embodiment explained using the above 
described flow chart in Fig. 27 has described the 
calculation of a data range whose maximum cone angle is 
narrow in the determining process of the operating data 
phase range by the operating data phase range 
determining means, but it is also possible to determine 
a data range for each voxel so that in step S4, the 
difference in the cone angle (corresponding to the 
absolute value of v) at the ends (data start/end 
positions) of the back projection data range becomes a 
minimum and perform a similar reconfiguration process 
based on this determined data range. 

An example of the method of determining the 
operating data phase range for each voxel (calculation 
of a data range whose back projection phase width is 
narrow) will be explained. First, the case where the z 
direction size (Z det ) of the radiation detector is 
sufficiently wide will be shown. When data can be 
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acquired in the same go-around phase range at all 
reconfigured voxels at the same z position (when 
reconfiguration is possible), or more specifically, 
when the image taking condition in Expression 51 is 
5 satisfied, the phase range having a small difference in 
the back projection phase range with respect to a voxel 
whose z position is located within the same plane is 
expressed by Expression 52, where 0 0 is the phase at" 
which the z position of the focus corresponds to the 
1.0 voxel position, dapp is the z direction size of the 
detector element and row is the number of detector 
rows . 

[Expression 51] 

j a p p » (row— 1) » (SOD-FOV/2) 

S I D r f , 0 . r FOV ^ 

— i7c + 2arcsin 

2 7c I [2 SOD JJ 

[Expression 52] 

0 o -f*/2^0<0 o +f7t/2 

However, when the relative moving speed 
between the examinee and the focus is high and it is 
15 impossible to acquire data at all voxels within the 
same phase range, that is, when the above described 
image taking condition is not satisfied, it is not 
possible to select the phase range as expressed in 
Expression 52. In such as case, it is possible to 
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determine the phase range using the method shown below 
If the distance between the radiation source 
and the rotation center is SOD, the relative movement 
distance of the radiation source relative to the 
examinee per rotation of the scanner on the radiation 
detector is J, the go-around phase of the fan beam 
source is p, the beam spreading angle between the beam 
directed to the reconfigured' voxel and central beam is 
a and the go-around phase of the parallel beam is (|), 
then the radiation source position S(p)=S(x s , y s , z s ) is 
expressed by Expression 2 described above. 

In terms of parallel beams obtained through , 
rearrangement process, this is expressed by Expression 
12 described above. 

Here, if the traveling direction of the 
parallel beam is w and the direction perpendicular to 
this w is t, then the t coordinate and w coordinate 
when the parallel beam with phase (j) passes through the 
coordinates (x, y) are expressed by Expression 13 and 
Expression 14 described above and the distance between 
the radiation source and tz plane (plane passing 
through the go-around axis and perpendicular to the 
parallel beam) is expressed by Expression 6 described 
above. Furthermore, when the parallel beam with phase 
passes through the reconfigured voxel (x, y, z) and 
crosses the detector whose distance from the radiation 
source is SID and the coordinates of the detector in 
the v axis (go-around axis) direction are H(x, y, ty) , 
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then this can be expressed by Expression 7 described 
above . 

Furthermore, if a phase range index is f, in 
order to back project the reconfigured voxel I (x, y, 
z) within a phase range having a small difference in 
the back projection phase range with respect to the 
voxel whose z position is located within the same 
plane, the z direction position of the radiation 
detector when the beams irradiated from the end 
positions B s and B e of the phase range fn used pass 
through the reconfigured voxel and cross the radiation 
detector must be located within the range of the 
detector, and therefore if the go-around phase when the 
z direction position of the focus is at the position of 
the reconfigured voxel is 9 0 , it is possible to select 
such (j) that satisfies Expression 53 and Expression 54 
and approximates to Q 0 -fn/2 infinitely. 
[Expression 53] 

H (x, y , <f>) ^ d a p p • (row— 1) / 2 
[Expression 54] 

H (x, y, f + f tc) ^— dapp (row— l)/2 

More specifically, when 9 0 =0, as in step S20 
shown in Fig. 31, if the initial value of <|> is -fn/2 and 
the phase accuracy to be calculated is Q (e.g., Q is a 
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phase angle by which the focus advances per one view) , 
if <|> is small ($<0) , H (x, y, (J>) decreases as <j> increases 
and H(x, y, <J>) increases as <J> decreases , and therefore 
if images are taken under reconf igurable conditions, 
the processes shown in Expression 55 and Expression 56 
are repeated until Expression 61 and Expression 62 are 
satisfied respectively as shown in steps S21 to S24. 
This makes it possible to satisfy Expression 57 and 
Expression 58 and select <|> so as to approximate to 0 0 - 
f7r/2 as much as possible. In this way, the phase range 
(Bs<(j)<Be) becomes the same as that expressed by 
Expression 12 and Expression 13. 
[Expression 55] 

if [ d a p p • (row— l)/2— H(x, y, <fr) < 0] , <*> = <f>+Q 
[Expression 56] 

if [dapp-(row— l)/2+H(x, y, <t> + f tz ) < 0] , <f> = <f> - Q 
[Expression 57] 

H(x, y, 0)^dapp - (row— l)/2 
[Expression 58] 



H(x, y, 0+f7c)^— dapp- (row— 1 ) / 2 



In this phase range calculation process, by 
determining the back projection phase range for each 
voxel so that the number of views is reduced, it is 
possible to improve time resolution for each voxel and 
obtain good image quality in regions where the examinee 
moves drastically by combining this with the 
aforementioned weighting back projection. Furthermore, 
by setting the back projection phase range for each 
voxel to within a time range in which images are taken 
at the same time wherever possible so that the time 
positions of the respective voxels in the displayed 
images come closer to one another, it is possible to 
shorten the time width contributing to the reconfigured 
image and improve time resolution. The back projection 
phase range in this case is ideally the same back 
projection phase range at all voxels, but even when the 
relative moving speed between the examinee and focus is 
high and it is impossible to obtain data at all voxels 
within the same phase range, it is possible to 
determine the back projection phase range for each 
voxel so that the examinee and focus come as close as 
possible to each other. 

To arbitrarily make changeable the 
relationship between a noise level and body axis 
resolution in the reconfigured image, a body axis (go- 
around axis) direction filter whose spatial frequency 
characteristic is changeable in the row direction is 
preferably superimposed on the projection data. This 
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superimposition of the body axis direction filter (body 
axis direction filtering) may be performed before or 
after the one-dimensional rearrangement process. The 
superimposition may also be included in the filter 
correction processing. Furthermore, the above 
described embodiment uses a tomograph using X-rays, but 
the present invention is not limited to such a 
tomograph and is also applicable to a tomograph using 
neutron rays, positron, gamma rays or light. 
Furthermore, the scan system is not limited to any one 
of the first-generation, second-generation, third- 
generation or fourth-generation systems, but can also 
be used for a multi-tube CT provided with a plurality 
of radiation sources, cathode scan CT or electron beam 
CT. Furthermore, the shape of the radiation detector 
is also applicable to any one of radiation detectors 
such as radiation detectors arranged on a cylindrical 
surface centered on a radiation source, plane 
detectors, radiation detectors arranged on a spherical 
surface centered on the radiation source, radiation 
detectors arranged on a cylindrical surface centered on 
the go-around axis, etc. Furthermore, the tomograph is 
not limited to a spiral orbit scan, but is also 
applicable to a circular orbit scan. Furthermore, 
projection data and reconfigured image that should 
actually be handled discretely are handled as 
continuous data, and therefore it is desirable to 
calculate discretely through interpolation in three 



directions of phase direction, row direction and 
channel direction of the radiation detector using an 
interpolation method such as Lagrange interpolation. 
Furthermore, the above described embodiment 
approximates arcsin with one approximate straight line, 
but it is also possible to approximate arcsin using a 
plurality of approximate straight lines (using 
different approximate straight lines according to the 
distance from the go-around axis) . Furthermore, a 
nonlinear function value of the present invention can 
also be calculated using advance calculations 
(tabulation) and interpolation for speed enhancement. 

In the above described embodiments, the 
process (S4) of determining the phase range of the 
projection data used for each voxel in Fig. 9 is also 
applicable to other embodiments of the reconfiguration 
means 22. 

The redundancy correction weighting process 
(S9) in Fig. 16 is also applicable to other embodiments 
of the reconfiguration means 22. 

The process (Sll) of rearrangement into data 
with a view of a multiple of 4 and process (S12) of 
grouping in Fig. 23 are also applicable to other 
embodiments of the reconfiguration means 22. 

The process (S13) of rearrangement in Fig. 24 
is also applicable to other embodiments of the 
reconfiguration means 22. 

The approximate straight line calculation 



process (S14) and parallel beam three-dimensional back 
projection process (S15) in Fig. 27 are also applicable 
to other embodiments of the reconfiguration means 22. 

As described above, according to the 
tomograph of the present invention, when 

reconfiguration is performed from data obtained through 
a scan, it is possible to reduce the distortion due to 
data discontinuity to a minimum and obtain images of 
high quality without producing any streak artifact in 
the reconfigured image. 

Furthermore, according to the tomograph of 
the present invention, it is possible to simplify an 
arcsin calculation used so far, enhance the speed 
drastically and obtain images of high quality in a 
short time by calculating an approximate straight line 
for a curve indicating the radiation source position 
with respect to the channel direction position of 
parallel beam projection data obtained by a parallel 
beam of a parallel shape viewed from the go-around axis 
direction generated from the radiation source. 

All the foregoing descriptions have been 
presented about the embodiments, but it is obvious for 
those skilled in the art that the present invention is 
not limited to these embodiments but can be altered or 
modified in various ways without departing from the 
spirit and accompanying claims. 

This application with claims of priority is 
based on Japanese Patent Application No . 2002-304 4 63 and 
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Japanese Patent Application No . 2003-078125, entire 
content of which is expressly incorporated by reference 
herein . 



